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ABSTRACT 


The  analysis  and  interpretation  of  the  reverberation  data  will 
require  attention  to  the  envelope  properties  of  the  narrowband 
process.  This  report  contains  preliminary  work  required  for  such 
a study.  The  work  was  performed  and  reported  at  this  time  due  to 
the  delay  in  data  collection. 


The  statistical  testing  of  the  sonar  return  data,  consisting  of 
individual  stave  outputs  from  several  hydrophone  arrays,  has  been 
delayed  due  to  the  inavailability  of  the  data  gathering  equipment, 
notably  the  computer  link  between  the  CDC  3200  and  the  Honeywell 
DDP-516.  This  link  has  been  in  heavy  use  during  the  past  quarter 
by  the  Signal  Physics  Division  which  is  involved  in  the  automatic 
classifier  program.  This  equipment  will  become  available  during 
the  next  quarter  and  the  stave  data  will  be  obtained  in  digital 
form. 
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To  detect  and/or  classify  target  echoes,  the  receiving  systems 
of  active  sonars  frequently  extract  and  process,  in  some  manner, 
the  envelope  of  the  received  signal,  which  may  be  reverberation  or 
reverberation  plus  the  echo  from  a target.  The  analysis  of  the  per- 
formance of  these  processors  requires  a minimum  amount  of  information 
about  the  statistical  descriptions  of  the  reverberation  envelope. 

To  a considerable  extent,  the  properties  of  the  envelope  are  determined 
by  the  characteristics  of  the  transmitted  waveform.  Thus,  the  covariance 
of  the  reverberation  process  is  very  similar  to  the  autocorrelation 
function  of  the  transmitted  signal.  The  signals  most  commonly  transmitted 
are  the  pulsed  sinusoid  (cv)  and  the  linearly  frequency  modulated  (FM) 
pulse.  Determining  the  properties  of  the  envelope  of  the  reverberation 
processes  generated  by  these  two  signal  types  is  the  goal  of  this  study. 

In  particular,  we  wish  to  know  the  covariances  and  intensity  spectra 
of  the  fluctuating  component  of  the  envelope , since  knowing  these  two 
second-order  statistics  is  necessary  for  the  analysis  of  a receiving 
system  that  makes  decisions  based  on  the  extracted  envelope. 

This  portion  of  the  study  is  organized  as  follows.  In  Chapter  II 
the  narrowband  reverberation  process  is  discussed,  with  emphasis  placed 
on  the  covariance  and  intensity  spectrum.  The  material  presented  here 
is  not  new,  but  it  is  needed  as  an  introduction  and  will  be  used,  in 
Chapters  III  and  IV,  in  the  computation  of  the  covariance  and  intensity 
spectrum  of  the  fluctuating  component  of  the  envelope  of  the  reverberation. 
In  Chapter  V,  the  results  obtained  thus  far  are  summarized  and  future 
research  in  this  direction  will  be  discussed. 
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II.  THE  NARROWBAND  REVERBERATION  PROCESS 


We  begin  by  assuming  that  the  reverberation  process  X(t),  defined 
on  the  time  (range)  interval  I = (T^  T2),  is  narrowband  (since  the 
transmitted  signals  that  generate  X(t)  are  assummed  to  be  narrowband) 
and  Gaussianly  distributed  with  zero  means.  Thus , 

(X(t))  = 0 , (2.1) 

where  <•>  denotes  the  infinite  ensemble  average.  Furthermore,  we 
assume  that  X(t)  is  a stationary  process  over  I,  an  assumption  generally 
not  true  in  practical  situations.  However,  for  reasons  of  analytical 
expedience,  this  assumption  is  made. 

The  narrowband  process  X(t)  can  be  written  as 


X(t)  = E(t)cos 


2irf  t + 
o 


(2.2) 


where  f , <J> ( t ) , and  E(t)  are,  respectively,  the  center  frequency,  the 
o 

phase,  and  the  envelope  of  X(t).  Since  X(t)  is  narrowband,  E(t)  and 
cos  4>(t)  are  slowly  varying  functions  when  compared  to  cos  2-irf^t. 

By  expanding  Eq.  (2.2),  X(t)  can  be  expressed  as 


X(t)  = x (t)cos  2?rf  t + x (t)sin  2mf  t , (2.3) 

c os  o 


where 


and 


x (t)  = -E(t)sin  $(t) 
s 

are  the  cosine  and  sine  quadrature  components  of  X(t). 


(2.1».b) 


The  covariance  of  X(t)  is  - , 

Kx(ti,t2)  i <X(t1)X(t2)>  • (2.5) 

When  X(t)  is  stationary,  then  K^t^tg)  = T = “ ^1* 

The  mean  intensity  (variance)  of  X(t)  is 

V = K^t.^)  = (x2(t1))  = ^(0)  » (2.6) 

which  is  independent  of  t under  the  assumption  of  stationarity . 


If  the  quadrature  representation  of  X(t)  (Eqs.  (2.U.a)  and  (2.U.b)) 
is  used  in  Eq.  (2.5),  then  the  covariance  can  be  written  as  (Ref.  1,  Sec. 9) 


K^t)  = 4'|po(t)cos  27rfQT  + Xq( i)sin  2TrfQTj  , 

where 

*P0(t)  S 2 (xc(tl)xe(t2))  + 2 (xs(tl)xs(t2)) 

and 

(VWV)  + 2 (‘.WV) 


(2.7) 


(2. 8. a) 


(2.8.b) 


b 


The  covariance  can  also  be  written  as 

K^t)  = Kq(t)cos  f 2iTfoT  + <t>o( x) j , (2.9) 

where 

Ko(T)  = Y[po2(t)  + Xo2(t)|1/2  * (2*10) 

and 


$o(t)  = tan-1  | Xo(t)/pq(t)]  , (2.11) 

are  the  envelope  and  phase,  respectively,  of  KX(T).  From  Eqs.  (2.8.b) 
and  (2.11)  we  see  that  <t>(o)  = 0,  and,  therefore,  using  Eq.  (2.9)  and 
Eq.  (2.6), 


III.  THE  FIRST  AND  SECOND  MOMENTS  OF  THE 
ENVELOPE  OF  THE  REVERBERATION 

This  study  is  primarily  interested  in  the  statistical  description 
of  the  envelope  E(t)  (see  Eq.  (2.2))  of  the  reverberation  process,  X(t), 
whose  partial  description  was  outlined  in  Chapter  II.  The  first  and 
second  moments  of  E(t)  are  (Ref.  2,  Sec.  9.1) 

/ \i/2 

<E(t)>  = ( \ y)  , (3.1) 

and 


^E(t)2^  = 2Y  , (3.2) 

where  V = (x(t)2)  is  the  variance  of  X(t).  Of  interest  here  is  the 
fluctuating  component  of  E(t) , which  is  the  variation  of  E(t)  about  its 
mean,  ^E(t))  , defined  by 

y(t)  = E(t)  - (E(t)>  , (3.3) 

with  (y(t))  = 0.  The  variance  of  y(t)  is,  using  Eqs.  (3.1)  and  (3.2), 

’ (r(tO 

= ^E2(t))  - (E(t))2 

m(2-f)v  ■ (3.U) 
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The  covariance  of  y(t)  is 


VW  - (y(ti)y<t2)> 


which,  using  Eq.  (3.3),  is 


(3.5) 


yw  - - <e<v>  <*y> 

-yw-f'  • 


(3.6) 


where 


WV  E (E(t1)E(t2)^ 


is  the  covariance  of  E(t).  If  we  assume  that  X(t)  is  a normal  random 
process,  then  we  can  use  a relationship  previously  obtained  for 
KE(t1,t2)  (Ref.  2,  Sec.  9-l), 


yw = I1 2y-i/2,-i/2;i,y) 


(3.7) 


where  F is  a hypergeometric  function  with  the  argument  k (t.  ,t  ) = k (t) 
f a i.  o x 2 o 

(see  Eq.  (2.13)).#  A general  series  form  of  is 

F (a  8'Y'x)  = 1 + — — 2L_  + a( ct+l) g( B+l)  2 ( ■*  p,\ 

2^!  ’P,Y’  ; 1 y II  y(y+1  21 '•*  • '3-8' 


Assuming  that  the  first  two  terms  of  this  series  accurately  represent  2F^ 
we  find,  with  Eqs.  (3.6),  (3.7),  (3.8)  that 


V W = Ky( t)  = J yk02('r)  ’ T = l2  " tl 


(3.9) 


•A  key  assumption  used  in  the  derivation  of  Eq.  (3-7)  is  that  the  envelope 
E(t)  has  a Rayleigh  first-order  probability  density  function;  we  also 
make  this  assumption. 


that  is,  the  covariance  function  of  the  fluctuating  component  of  the 
envelope  of  the  reverberation  is  proportional  to  the  square  of  the 
envelope  of  the  covariance  function  of  the  reverberation  process  X(t). 
This  result,  with  comparisons  to  experimental  data,  is  given  by 
Ol'shevskii  in  Ref.  3,  Ch.  IV. 

The  expression  for  K (t),  Eq.  (3.9),  is  approximate  since  only 
the  first  two  terms  of  the  series  of  Eq.  (3.8)  were  retained.  To 
determine  the  amount  of  energy  lost  in  this  approximation  we  note 
that  the  time  variance  of  y(t)  is  given  by  Eq.  (3.M, 


ay2  = (2  - lOY  * 0.43  V 


2 

whereas  our  approximation  of  a is  (Eq.  ( 3 . 9 ) ) 


vo> 


approx 


= -g4'  = 0.394'  = 0.91  a 2 

o y 


Therefore,  about  10%  of  the  energy 
following  the  second  of  the  series 


of  y(t)  is  contained  in  the  terms 
of  Eq.  (3.8). 
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IV.  THE  INTENSITY  SPECTRUM  OF  THE  FLUCTUATING 
COMPONENT  OF  THE  ENVELOPE 

Applying  the  Wiener-Khint chine  (W-K)  theorem  (Ref.  2,  p.  lUU)  , 
the  intensity  spectrum  of  the  fluctuating  component  (see  Eq.  (3.3)) 
can  be  obtained  from  the  covariance,  K^(t),  of  y(t): 

W (f)  = U / K (t)cos  urdT,  w = 2irf  . (U.l) 

y JQ  y 

Using  the  approximation  Eq.  (3.9)  for  K (x),  the  intensity  spectrum 
of  y(t)  is 


U/-»r  / O 

Wy(f)  = k^  (x)cos  wTdT  . (U.2) 

The  explicit  form  of  k^f)  is  determined  mainly  by  the  characteristics 
of  the  waveforms  of  the  transmitted  signals  that  generate  the  reverbera- 
tion process  X(t).  The  covariance  ■'■s  approximately  equivalent 

to  the  autocorrelation  function  of  the  transmitted  signal  (Ref.  h,  Ch.  V) 


s(t)s(t+i)at 


(h.3) 


where  S(t)  is  the  narrowband,  transmitted  signal.  Under  certain 
conditions,  such  as  no  scatterer  Doppler  or  platform  motion,  local 
stationarity,  etc. , we  can  write  (within  a scale  factor) 

K^t)  = C(t)  . (it.1*) 


f. 
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Before  proceeding,  the  transmitted  signal  must  be  specified.  Currently 
the  two  types  of  signals  that  are  most  often  transmitted  by  an  active 
sonar  are  the  cw  pulse  and  the  linearly  frequency  modulated  (FM)  waveform 
Each  of  these  will  be  considered  in  the  next  progress  report. 


V.  CONCLUDING  REMARKS 


With  the  necessary  first-  and  second-order  statistics  of  the 
reverberation  envelope  we  can  go  further  and  determine  the  properties 
of  the  reverberation  envelope  that  has  been  passed  through  various 
types  of  linear  systems.  One  problem  to  be  solved  is  that  of 
determining  the  statistical  properties  of  the  envelope  that  has  been 
subjected  to  a smoothing  process,  i.e.,  averaged  with  a fixed-length, 
lagged  time  window.  Of  considerable  interest  here  is  the  comparison 
of  the  FM  and  cw  envelope  processes  that  have  been  subjected  to  the 
same  linear  system.  Since  the  input  spectra  of  these  two  processes 
are  different,  it  can  be  assumed  that  the  corresponding  output  spectra 
will  be  different.  In  the  case  of  smoothing,  a variable  parameter  is 
the  integration  time  which,  when  increased,  results  in  a decrease  in 
the  spectral  widths  of  the  envelope  processes,  and  hence  a decrease 
in  the  background  reverberation  level  in  the  sonar  receiver. 
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